AD-AQ74  957  RENSSELAER  POLYTECHNIC  INST  TROY  N  Y  DEPT  OF  MECHANI— ETC  F/G  20/11 

A  THEORY  OF  THERMOVISCOPLASTICITY  FOR  UNIAXIAL  MECHANICAL  AND  T-ETC(U) 
JUL  79  E  P  CERNOCKY*  E  KREMPL  N00014-76-C-0231 

UNCLASSIFIED  RPI-CS-79-3  NL 


I  OF  | 

*0^4  95  7 


l< 


fui^lfrwipniwi  [*iihjwhi»» 


A  THEORY  OF  TOERMOVXSCO PLASTICITY  FOR 
UNIAXIAL  MECHANICAL  AND  THERMAL  LOADING 

E.P.  Cemocky  and  E.  Krempl 
Department  of  Mechanical  Engineering, 
Aeronautical  Engineering  &  Mechanics 
Rensselaer  Polytechnic  Institute 
Troy,  New  York  12181 


Report  No.  RPI  CS  79-3 
JUly  1979 


has  boon  Gpprovod 

idf  public  rslsase  and  im 
distribution  is  unlimited. 


A  THEORY  OP  THERMDVISCO PLASTICITY  FOR  UNIAXIAL 
MECHANICAL  AND  THERMAL  LOADING 


E.P.  Cemocky  and  E.  Krempl 
Department  of  Mechanical  Engineering, 

Aeronautical  Engineering  &  Mechanics 
Rensselaer  Polytechnic  Institute 
Troy,  New  York  12181 

ABSTRACT 

A  previously  proposed  three-dimensional  isotropic  theory  of  thermo¬ 
viscoplasticity  based  on  total  infinitesimal  strain  is  specialized  to  a  uni¬ 
axial  state  of  stress.  This  uniaxial  theory  consists  of  a  first-order  dif¬ 
ferential  constitutive  equation  linear  in  the  mechanical  strain  rate  and  the 
stress  rate  but  nonlinear  in  the  mechanical  strain,  the  temperature,  and  the 
stress.  This  equation  is  coupled  to  a  constitutive  heat  equation  where  the 
work  in  a  homogeneous,  adiabatic  cycle  which  starts  and  ends  at  zero  stress 
is  completely  converted  into  temperature  change.  In  cyclic  loading  the 
mechanical  constitutive  equation  is  augmented  through  a  procedure  which  we 
call  "storage  and  updating". 

The  qualitative  solution  properties  of  this  system  of  differential 
equations  are  investigated.  Assuming  adiabatic  conditions  it  is  shown  that 
near  the  stress-strain- temperature  origin  or  under  large  instantaneous 
changes  in  the  strain  (stress)  rate,  the  predicted  material  behavior  is 
thermoelastic.  The  solutions  for  large  time  under  monotonic  loadings  are 
given.  A  set  of  material  constants  is  assumed;  the  equations  are  numerically 
integrated  to  simulate  homogeneous  monotonic  loading  and  cycling  and  corre¬ 
sponding  deformation-induced  temperature  changes. 

Other  examples  include  thermal  monotonic  and  cyclic  loading  under 
constraint  (thermal  fatigue  cycling)  and  mechanical  cycling  with  simultaneous 
rapid  heating. 


Introduction 


Mechanical  deformation  is  often  idealized  to  occur  under  isothermal 
conditions.  This  idealization  is  a  convenient  and  a  reasonable  simplifying 
assumption  under  many  conditions.  In  reality  mechanical  deformation  induces 
temperature  changes  in  materials,  under  normal  (slow)  rates  of  loading 
these  deformation-induced  temperature  changes  go  unnoticed  due  to  the  equi¬ 
librating  effect  of  heat  conduction.  However  under  conditions  of  rapid 
monotonic  loading  or  under  rapid  cyclic  loading  at  sufficiently  high  stress 
levels,  inelastic  self-heating  of  materials  can  become  important  and  may 
result  in  melting  of  the  material  [1] .  In  particular  situations  it  may 
become  necessary  to  calculate  accurately  the  temperature  changes  generated 
during  the  course  of  mechanical  deformation;  examples  are  given  in  [2]. 

Here  we  present  a  coupled  theory  of  thermoviscoplasticity  for  the  uni¬ 
axial  state  of  stress  and  for  infinitesimal  deformations.  To  represent  the 
deformation-induced  temperature  changes  which  occur  during  mechanical  de¬ 
formation,  we  postulate  a  constitutive  heat  equation.  Hiis  equation  is 
coupled  to  a  differential  equation  governing  the  mechanical  deformation 
Which  is  called  the  mechanical  constitutive  equation.  Both  equations 
share  a  common  nonlinear  dependence  upon  the  present  stress,  strain,  and 
temperature,  and  employ  total  strain  only  without  any  decomposition  of 
strain  into  separate  elastic  and  inelastic  components. 

Recent  experimental  studies  [3-6]  have  indicated  that  significant 
rate  dependence  can  be  present  in  the  mechanical  behavior  of  metals  at 
elevated  and  at  room  temperatures.  Ihe  mechanical  constitutive  equation 
proposed  herein  is  intended  to  represent  both  this  rate-dependent  behavior, 
and  with  augmentations  the  hysteretic  behavior  observed  in  metals,  including 
cyclic  hardening,  cyclic  softening,  and  the  Bauschinger  effect.  Under  the 
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idealization  of  isothermal  deformation  the  constitutive  equation  reduces  to 
the  previously  proposed  isothermal  mechanical  constitutive  equation  [7] . 

The  purpose  of  this  paper  is  to  examine  the  behavior  predicted  by  the 
proposed  constitutive  equations  under  homogeneous  deformations  using  a 
variety  of  different  loading  histories  including  mechanical  and  thermal 
cycling.  This  is  accomplished  by  examining  the  behavior  of  the  proposed 
nonlinear  coupled  differential  equations  first  by  analytical  means.  For 
hypothetical  but  realistic  material  properties  the  coupled  equations  are  then 
numerically  integrated  using  the  Bulrisch-Stoer  Algorithm  [22]  to  illustrate 
the  solution  properties  in  stress-strain  and  temperature-strain  (stress)  graphs. 

The  Constitutive  Equations 

We  designate  o  and  e  as  the  axial  stress  and  infinitesimal  strain, 
respectively;  6  is  the  absolute  temperature  and  0q  is  the  equilibrium  or 
reference  temperature.  Young's  modulus  is  E[0]  where  throughout  this  paper 
brackets  denote  function  dependence  upon  the  indicated  argument.  The  coef¬ 
ficient  of  thermal  expansion  is  denoted  as  a  and  the  specific  heat  of  the 
material  is  C;  here  a  and  C  may  depend  upon  absolute  temperature  or  may  be 
taken  as  constants.  A  superposed  dot  denotes  time  differentiation,  while 
the  usual  comma-notation  is  used  to  represent  differentiation  with  respect 
to  position  coordinates.  The  constant  density  of  the  material-  is  denoted 
as  p. 

We  introduce  the  uniaxial  mechanical  strain  Z  as 

Z  ■  e  -  ar(0  -  0  )  (1) 

o 

Following  [6]  we  propose  our  uniaxial  mechanical  constitutive  equation  in 


its  simplest  form 
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E[0]z-  a-  (a-g[Z,0))/k[CT-g[z,e]J  .  (2) 

The  function  k[  ]  is  a  continuous,  positive,  even,  and  decreasing  function 
of  its  argument.  The  function  g[Z,  0]  represents  the  equilibrium  stress- 
strain  curve  of  the  material  at  temperature  0.  As  defined  in  [7-10]  the 
function  g[  ]  is  odd  in  Z  and  for  fixed  0  and  positive  Z  g[  ]  has  the 
appearance  of  a  tensile  stress-strain  curve;  g  is  initially  linear  elastic 
in  Z.  Hie  second  argument  of  g[  ]  represents  the  temperature  dependence 
of  material  parameters  contained  in  g [  ]  such  as  E [ 0] .  He  denote  the  final 
tangent  modulus  (in  Z)  of  g[  ]  as  £^[0],  so  that  for  a  sufficiently  large  Z 
we  have  g'  =  ^  ^  «  E^IO]*.  The  differential  equation  (2)  is  linear  in 

the  stress  rate,  strain  rate,  and  temperature  rate,  but  it  is  nonlinear  in 
the  stress,  strain,  and  temperature  through  the  nonlinear  functions  k[  ] 
and  g[  ]. 

Within  the  context  of  infinitesimal  theory  the  first  law  of  thermo¬ 
dynamics  may  be  written  as 

pe  »  ae  -  qA  A  +  pR  (3) 

where  e  is  the  internal  energy  per  unit  mass,  q.  .is  the  divergence  of 
the  heat  flux  vector,  and  R  is  an  internal  heat  supply  (which  is  externally 
specified  and  controlled;  for  example,  electric  resistance  heating).  We 
postulate  that  the  time  rate  of  change  of  the  internal  energy  is  equal  to 
the  following  exact  differential  in  stress  and  temperature  . 


* 

Here  we  assume  the  stress-strain  diagram  is  ultimately  linear  in  Z.  This 
assumption  is  convenient  but  not  necessary  for  the  theory. 


pe 


pC[0]9  +  dt  (2Er0] 


+  a[0]ae) 


(4) 


Combining  the  constitutive  assumption  (4)  with  the  first  law  of  thermo¬ 
dynamics  (3)  gives  us  our  equation  of  heat  conduction* 

(pC  +  oa)0  +  q.  .  -  ae  -  (a  +  Ear0)  ^  +  pR  (5) 

1^1  E 

where  in  (5)  we  surpress  writing  the  temperature  dependence  of  E[  ],  C[  ], 
and  art  1 .  To  examine  heat  conduction  one  must  additionally  postulate  a  con¬ 
stitutive  choice  for  the  heat  flux  vector  q.  Since  we  will  consider  only 
homogeneous,  adiabatic  deformations  in  examining  the  predictions  of  this 
theory,  q^  ^ =  0  and  a  specific  choice  for  the  heat  flux  vector  is  not 
necessary. 

We  note  that  the  constitutive  assumption  (4)  allows  us  to  write  in 
adiabatic  deformation  without  heat  supply 

e-e[o[t],0[t],0o].  (6) 

For  adiabatic  loading  and  for  zero  internal  heat  supply,  Eq. (3) 
indicates  that  all  of  the  mechanical  work  must  be  converted  into 
change  of  the  internal  energy.  Experiments  indicate  that  this  energy  change 
is  manifested  through  a  small  percentage  of  the  mechanical  work  («  10%) 
being  converted  into  microstructural  change  and  the  remaining  major  portion 
of  the  mechanical  work  being  converted  into  heat  [11-14],  For  the  purposes 
of  this  paper  we  assume  that  the  mechanical  work  can  be  completely  transformed 
into  temperature  change;  this  assumption  corresponds  to  Eq. (5) 


* 

Here  and  throughout  this  paper  we  assume  that  quantities  containing  time 
derivatives  of  temperature-dependent  material  parameters  can  be  neglected. 
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To  examine  the  predictive  capability  of  the  theory  we  first  assume  the 
absence  of  any  externally  imposed  internal  heat  supply  (R«0)  in  adiabatic 
deformation.  Ihe  heat  equation  (5)  may  then  be  rewritten  as 


(pC+OCT)  *CT-  (a+BQf0)|f/E 


or  equivalently  using  (2)  we  have 

**  - « i£^rtiU  - 


or 


Pce  -  .  Eaefe 


(7) 


(8) 


(9) 


In  all  of  Eqs.  (2)  -  (9)  we  may  choose  either  E«E[6]  or  E«  constant. 


Cyclic  Loading 

In  monotonic  tensile  or  compressive  deformations,  the  sign  of  the 
overstress  (o-g)  predicted  by  Eqs.  (2)  and  (7)  does  not  change  (non-negative 
in  tension,  nonpositive  in  compression) .  However  in  cyclic  deformations 
the  predicted  stress-strain  response  curve  can  be  forced  to  cross  the 
equilibrium  curve  g[  ]  during  unloading  and  reloading*.  When  such  crossing 
occurs  the  overstress  (<j-  g)  vanishes  and  subsequently  changes  sign.  At 
points  where  the  overstress  vanishes  we  (discontinuously)  augment  the  equi¬ 
librium  function  g[  ]  in  order  to  recognize  the  microstructural  changes  that 
have  occurred  in  deforming  the  material  up  to  this  point.  Augmentation 
consists  of  a  process  which  we  call  storage  and  updating. 
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This  cannot  occur  due  to  creep  or  relaxation 


At  the  occurrence  of  vanishing  overstress,  we  store  the  values  of  the 
stress,  strain,  and  temperature.  We  simultaneously  may  change  the  value  of 
material  constants  in  the  function  g[  ]  to  complete  our  representation  of 
prior  micros true tural  change.  The  constant  Xf  is  most  suitable  for  this  and 
in  this  paper  the  value  of  Xf  is  updated  to  a  new,  experimentally  determined 
value*.  We  then  continue  the  deformation  using  the  augmented  equilibrium 
function  in  the  constitutive  equations. 

The  general  idea  about  augmentation,  i.e.,  the  representation  of  the 
microstructural  changes, was  given  in  [15]  and  the  manner  and  degree  of  change 
of  parameters  has  been  presented  for  the  idealization  of  rate- independent 
materials  in  [16].  The  procedure  used  in  [16]  for  the  determination  of 
parameter  evolution  can  be  extended  to  this  case  of  rate-dependent  materials 
if  the  stress-strain  curve  of  [16]  is  identified  as  the  equilibrium  stress- 
strain  curve  g[  ]  of  this  report. 

Upon  initial  deformation  the  material  is  represented  through  the  g- 
function  which  contains  material  parameters  such  as  Xf  [10],  which  has  a  con¬ 
stant  initial  value  denoted  as  xj0'.  We  express  this  state  of  the  material 
by  writing 

g-g[z, 8]  |  (Q)  do) 

xf 

As  the  material  is  cycled  the  overstress  vanishes  prior  to  changing 
sign.  At  this  instant  we  denote  the  values  of  the  stress,  strain,  tempera¬ 
ture,  and  mechanical  strain  as  and  Z^,  respectively, 

and  we  store  these  constants.  We  simultaneously  update  the  parameter  Xf 

* 

The  parameter  X^  here  refers  to  the  g-function  representation  in  [10] .  if 
necessary  additional  material  parameters  can  be  updated. 
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(D* 

by  changing  it  to  a  new  value  Xf  ,  and  we  obtain  our  augmented  equilibrium 
function  by  replacing  the  representation  of  Eq. (10)  with 

g  =  g[Z- Z(1),0]  +  a(1)  (11) 

v '  1 

xf 

We  continue  the  deformation  using  (11)  until  the  next  occurrence  of  vanishing 
overstress,  where  the  process  is  repeated.  Hence  at  the  nth  occurrence  of 
vanishing  overstress  we  obtain  through  the  storage  and  updating  procedure 

g  =  g[Z-  Z(n),6]  (r)  +  a(n)  (12) 

Xf 

To  represent  the  case  of  cyclic  hardening  the  parameter  Xf  increases 
While  in  cyclic  softening  Xf  decreases.  In  the  case  of  steady-state  cyclic 
response,  the  storage  process  continues  but  the  parameter  Xf  ceases  to 
change  and  remains  at  a  fixed  value. 

Through  this  storage -updating  procedure  the  material  "remembers"  the 
values  of  the  stress,  strain,  and  temperature  at  each  occurrence  of  vanish¬ 
ing  overstress.  Between  occurrences  of  vanishing  overstress  the  micro- 
structural  changes  are  incorporated  into  the  model  through  the  nonlinear 
character  of  g[  ]  and  of  k[  ]  and  through  the  storage  of  the  stress,  strain, 
and  temperature  at  all  points  of  vanishing  overstress.  Although  this  storage¬ 
updating  procedure  is  applied  discontinuously,  the  response  predicted  by  the 
model  remains  smooth.  This  procedure  is  independent  of  the  particular  control 
used,  i.e.,  whether  stress,  strain,  or  temperature  is  prescribed.  The  details 
of  the  augmentation  procedure  will  be  described  in  a  forthcoming  paper. 

* 


The  new  value  must  be  known  from  experiments,  see  Pig. 7  of  [16] 


Figure  1  shows  the  hyrfc  uresis  loop  obtained  through  the  numerical  solu¬ 
tion  of  the  coupled  equations  of  this  model  using  this  storage-updating 
procedure;  the  input  is  a  1  Hz  strain-controlled  sine  wave,  and  the  functions 
and  material  parameters  used  appear  in  Table  1.  Figure  2  shows  the  predicted 
temperature  response  corresponding  to  the  adiabatic  cycling  of  Figure  1. 

Both  elastic  cooling  and  inelastic  heating  occur  over  each  cycle  with  a 
resultant  net  increase  of  temperature  with  each  cycle.  The  qualitative 
features  of  Figure  2  correspond  to  experimentally  observed  metal 
behavior  [14] . 


Behavior  in  the  Neighborhood  of  the  Stress-Strain-Temperature  Origin 

To  model  metal  deformation  the  function  g[Z, 0]  should  be  approximate1 y 
linear  thermoelastic  so  that 


and 


g[Z,  0]  «  E[0]  Z  when  Z  «  0 
E[0]  E[0  ]  when  0^0 


(13) 

(14) 


The  condition  (13)  may  be  achieved  through  the  representation  of  the  func¬ 
tion  g[  ]  proposed  in  [10],  When  (13)  applies  the  differential  equation  (2) 
reduces  to 

a  «  EZ  (15) 

and  v  e  initial  mechanical  behavior  is  linear  thermoelastic. 

Similarly  the  heat  equation  at  c =  0  reduces  to 


d0 

de 


-EcyQ 


pc  -  Ecr  9 


-  Edf0 
pc 


(16) 


and  the  initial  temperature  response  is  thermoelastic  [in  Eq. (16)  the  term 
2 

Ecr  0  is  negligible  relative  to  pC] .  Equation  (16)  indicates  that  the  initial 


temperature  response  is  independent  of  the  value  of  the  strain  rate  or 
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stress  rate  and  is  cooling  in  tension  and  heating  in  compression.  Subse¬ 
quently  thermoinelastic  heating  follows  when  the  mechanical  behavior  becomes 
inelastic. 

We  note  that  the  abo,Te  properties  apply  also  in  the  neighborhood  of 
vanishing  overstress. 

Limiting  Behavior  Predicted  by  the  Constitutive  Theory 

To  examine  the  behavior  predicted  by  this  constitutive  theory  in 
constant  strain  rate  and  constant  stress  rate  deformations  at  the  extremes 
of  loading  rates  and  at  large  times,  we  rewrite  the  differential  equation  (2) 
as  a  nonlinear  integral  equation,  see  [8] .  Using  the  methods  of  [7-9]  we 
perform  the  limit  of  the  stress  response  predicted  under  limitingly  fast 
and  slow  constant  strain  rate  deformation  and  obtain 


and 


lim  g  =  g[Z,  0] 
e-»  0 


lim 

e  — •  oo 


dg 

de 


=  E[0] 


(17) 


(18) 


Corresponding  to  (17)  and  (18)  the  heat  equation  (7)  reduces  to  the  re¬ 
spective  forms 


d9 

de 


£(■-(!— )*« 


(19) 


and 


d0  Eor0 

de  pC 


(20) 


where  in  (19) 


_d _ 3_  ,  d£  _a_ 

de  “  de  de  90 


(21) 


* 


For  constant  E  (18)  reduces  to  g=EZ. 


Equations  (17)  and  (18)  represent  the  limiting  lower  and  upper  bounds, 
respectively,  upon  the  mechanical  response  obtained  in  constant  strain-rate 
deformation;  the  same  limits  are  obtained  for  constant  stress-rate  loading 
[7,8],  It  is  because  of  the  limit  (17)  that  we  identify  the  function  g[Z, 0] 
as  the  equilibrium  response  of  the  material.  From  the  heat  equation  in  the 
case  of  infinitely-slow  loading  we  observe  that  the  temperature  response  is 
initially  thermoelastic  but  subsequently  is  thermoinelastic.  When  the  equi¬ 
librium  response  g[  ]  ceases  to  be  elastic  in  Z,  Eq. (19)  predicts  inelastic 
self-heating  in  the  material.  The  rate  of  temperature  change  in  this  in¬ 
elastic  heating  is  less  than  that  which  occurs  at  nontrivial  strain  (stress) 
rates.  In  the  case  of  infinitely  fast  loading,  the  mechanical  behavior  is 
thermoelastic,  and  the  corresponding  heat  equation  predicts  only  thermo¬ 
elastic  behavior  -  tensile  cooling  and  compressive  heating. 

In  a  constant  strain  rate  deformation  we  take  the  limit  t-*«o  to  deter¬ 
mine  the  limiting  or  steady-state  value  of  the  overstress  (a-  g)  and  we 
denote  this  limit  as 

{ex-  g}  =  lim  (a-  g)  (22) 

t~»oo 

Following  [7-9]  we  use  the  integral  representation  of  the  mechanical  consti¬ 
tutive  equation  to  determine  this  limit  and  we  obtain 

{a-g}=  (E-Efc)(l-a  ||)ek[{a-g}]  (23) 

Equation  (23)  is  a  transcendental  equation  in  {a-g}  which  indicates 
that  the  steady-state  overstress  depends  nonlinearly  upon  the  applied  strain 
rate.  Alternatively  in  constant  stress  rate  deformation  the  differential 
equation  (2)  is  substituted  into  the  integrand  of  the  integral  equation  fol¬ 
lowing  [8]  to  obtain  the  limiting  overstress  in  constant  stress  rate 


deformation 


£cr  -  g}  =  (E-  Et)k[{a-  g}]a/Efc  . 


(24) 


Hie  equivalence  of  the  limits  (23)  and  (24)  is  established  through 
execution  of  the  time  limit  of  the  time  derivative  of  the  integral  repre¬ 
sentation  which  furnishes 


lim 

t-*oo 


da 

de 


■*(*-•  S) 


(25) 


in  constant  strain  rate  or  constant  stress  rate  deformation. 

The  limits  (23)  -  (25)  are  rapidly  attained  asymptotic  limits  of  the 
response  predicted  by  the  constitutive  theory.  Although  we  let  t  -*  ®  to 
obtain  the  limits,  the  mechanical  response  rapidly  locks-into  the  steady- 
state  response  corresponding  to  these  limits  at  very  small  strains,  as  shown 
in  Pigs. 3  and  5  and  the  examples  in  [7-9], 

Numerical  studies  of  this  constitutive  theory  in  simulated  deformations 

d8 

indicate  that  the  term  a  —  appearing  in  (23)  is  negligible. 

ae 

Both  (23)  and  (24)  indicate  that  the  difference  between  steady-state 
flow  stress  and  steady-state  equilibrium  response  depends  nonlinearly  upon 
the  applied  strain  rate  or  stress  rate.  Because  of  this  fact,  in  deforma¬ 
tions  in  which  the  strain  (stress)  rate  is  increased  by  many  orders  of  mag¬ 
nitude  as  in  Figs. 3  and  5,  the  difference  between  stress  response  and  equi¬ 
librium  response  in  steady-state  behavior  may  not  increase  by  several  orders 
of  magnitude;  a  small  increase  in  this  steady-state  spacing  can  be  modeled 
if  an  appropriate  k[  ] -function  is  used.  Consequently  increasing  the  strain 
(stress)  rate  can  result  in  a  significant  but  small  increase  in  steady-state 
flow  stress  which  is  controlled  by  k[  ].  This  prediction  qualitatively  cor¬ 
responds  to  the  behavior  reported  in  experiments,  see  [3-5] . 


12 


I 

I 

I 

I 

I 


Behavior  Predicted  Under  Jumps  in  Strain  Rate  or  Stress  Rate 

We  examine  the  behavior  predicted  under  instantaneous  jump  increases 
(decreases)  in  the  applied  strain  (stress)  rate.  We  let  e  denote  a  tensile 
strain  rate  prior  to  a  jump  increase  or  decrease  in  the  strain  rate  and  we  let 
e+=Y®~  represent  the  strain  rate  jump;  e+  denotes  the  strain  rate  after  the 
jump  (we  similarly  use  +  and  -  in  the  case  of  the  stress-strain  and  temperature- 
strain  slopes).  Hie  factor  y  is  a  large  positive  constant  in  the  case  of  a 
strain  rate  jump  increase  and  a  small  positive  constant  in  the  case  of  a  strain 
rate  jump  decrease  in  which  loading  is  continued;  in  the  case  of  a  strain  rate 
jump  decrease  in  which  reversal  (unloading)  occurs,  y  is  negative.  From 
Eq. (2)  and  the  chain  rule  we  have 


=  E  _  CT-g[Z.9] 


de 


de  kf  ]t 


(26) 


Using  (26)  both  before  and  after  the  jump  in  strain  rate,  we  obtain  at 
the  point  of  jump  in  rate 

+  ._+ 

+ 

Y 


52_  =  E  + 1  (Aa~  .  p  *z~\ 

de  de  y  \  de  de  / 


(27) 


which  relates  the  stress-strain  slopes  before  and  after  the  jump  in  strain 
rate.  Proceeding  similarly  with  Eq. (7)  we  obtain 
,+ 


de  1  de  1  r/da  1  da”Y  /  1M 

de  "  y  de  (pC  +  aa)  L\ de  “  y  "de/ (a  +  Ea0)/E  _  a(1  ”  y/J  * 


(28) 


When  a  large  jump  increase  in  the  strain  rati*  is  imposed  y  »  1  (typically 
2  1 

Y  "  10  )  and  the  terms  containing  -  become  negligible  in  (27)  and  (28)  we  obtain 


da+  /.  d6  \ 

-X—  **  E(  1  -  a  -~)«E 
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We  observe  that  the  mechanical  and  thermal  responses  are  approximately  linear 
thermoelastic  in  slope  in  departing  the  stress-strain- temperature  point  cor¬ 
responding  to  strain  rate  jump  increase.  The  mechanical  response  upon  strain 
rate  jump  increase  is  shown  in  the  numerical  simulation  of  Pig. 3  While  the 
corresponding  temperature  response  is  shown  in  Pig. 4.  In  Fig. 4  we  observe 
thermoelastic  cooling  immediately  following  the  strain  rate  increase 
corresponding  to  the  prediction  of  (30). 

Alternatively  if  we  consider  a  jump  decrease  in  the  strain  rate  in  the 

region  of  steady-state  stress  response  «  E^  while  continuing 

-2  1 

tensile  loading,  then  1  »  y  >  0  (typically  y  =  10  )  and  the  terms  with  —  become 

dominant  in  (27)  and  (28).  We  obtain  the  prediction  of  a  very  large  negative 
slope  for  the  stress  response  departing  the  point  of  strain  rate  jump  decrease 
^here  E/y)  as  shown  in  Pig. 3*.  This  rapid  decay  of  stress  response 

results  in  strong  heating  in  the  material,  as  predicted  by  (28)  where  1/y 
dominates,  and  this  heating  is  illustrated  in  Fig. 4  corresponding  to  the 
strain  rate  jump  decrease  of  Fig. 3. 

In  the  cases  of  jump  increase  in  strain  rate  and  jump  decrease  in 
strain  rate,  the  stress  response  asymptotically  merges  with  the  response 
curve  which  would  have  been  obtained  under  initial  uniform  application  of 
the  high  (low)  constant  strain  rate  deformation;  see  Fig. 3.  Thus  this  theory 


We  emphasize  that  this  decrease  does  not  represent  a  reversal,  but  con¬ 
tinued  loading. 
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predicts  an  absence  of  "the  strain  rate  history  effect"  upon  stress  response; 
this  prediction  corresponds  to  the  experimental  behavior  reported  in  [3,4]  . 
However,  a  "strain  rate  history  effect"  upon  temperature  response  under  jump 
in  strain  rate  is  evident,  Fig. 4;  experiments  have  not  yet  been  performed  to 
test  for  this  effect. 

To  study  the  behavior  under  jump  in  stress  rate  we  use  the  chain  rule 
to  rewrite  (2)  as 

A  f*  1  *9  Q  rr  _  r» 

(31) 


de  _  1  .36  ,  cr-  g 
da  E  da  Ekf]a 
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Following  the  previous  example  we  let  a  =  0a  and  we  consider  jump 


increase  (Q  »  1)  and  jump  decrease  (0  «  1)  in  the  stress  rate  in  the 

da" 


steady-state  region  of  stress  response  where  «  E.  Using  (31)  both 
before  and  after  the  jump  in  stress  rate  we  obtain  at  the  point  of  jump  in 
stress  rate 


de*  1  de*  ,1  ( de"  _  1  _  d£\ 

da  E  da  0  V da  E  d o' 


(32) 


while  from  the  heat  equation  we  obtain 
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(33) 


Equations  (32)  and  (33),  respectively,  relate  the  strain-stress  and  the 


temperature-stress  slopes  before  and  after  the  jump  in  stress  rate.  With  a 


sufficiently  large  jump  increase  in  stress  rate  the  terms  involving  ^  are 


negligible,  e.g.,  if  ^  and  we  obtain 


In  the  case  of  a  hypothetical  material  extremely  sensitive  to  temperature 
changes,  a  strain  rate  history  effect  would  be  predicted  due  to  different 
temperatures  resulting  at  the  same  strains  in  the  different  loading  histories. 
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and  using  (34)  in  (33) 
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We  observe  that  under  a  large  instantaneous  jump  increase  in  stress  rate, 
both  the  mechanical  and  thermal  response  are  approximately  linear  thermo¬ 
elastic  in  departing  the  point  of  jump  in  stress  rate;  this  behavior  is 
illustrated  in  the  numerical  simulations  of  Figs. 5  and  6. 


Alternatively  if  we  consider  a  jump  decrease  in  stress  rate  then  ^  is 


_2 

very  large  (typically  0»10  )  and  from  (32)  we  obtain 

de  de  ‘ 


(36) 


Thus  under  large  jump  decrease  in  stress  rate  the  mechanical  response  does 

not  have  a  large  negative  slope;  rather  the  response  departs  the  point  of 

jump  decrease  in  stress  rate  with  a  very  shallow  stress-strain  slope,  as 

da+ 

shown  in  Fig. 5.  Indeed  — —  may  become  nearly  zero  for  very  large  jump 

de 

decreases  in  stress  rate  (very  small  positive  0).  This  indicates  a  strong 
bias  in  the  theory  between  behavior  under  strain  control  with  strain  rate 
jump  decrease  and  behavior  under  stress  control  with  stress  rate  jump  de¬ 
crease;  such  a  bias  has  been  observed  in  experimental  studies  [3,4].  In 
the  case  of  strain  rate  and  stress  rate  jump  increase  the  theory  does  not 
predict  any  bias. 

From  (33)  in  the  case  of  jump  decrease  in  stress  rate  we  again  obtain 
a  prediction  of  large  rapid  heating  due  to  the  dominance  of  i  in  (33)  and 
this  prediction  is  illustrated  in  the  numerical  simulation  of  Fig. 6  which 
corresponds  to  the  mechanical  behavior  of  Fig. 5. 


Thermal  Monotonic  Loading  and  Thermal  Cycling 

The  theory  is  capable  of  representing  loadings  in  which  both  the  me¬ 
chanical  quantities  (stress  or  strain)  and  the  temperature  are  controlled. 

This  requires  specification  of  the  internal  heat  supply  term  R(t)  and  here 
we  consider  large  R(t)  in  order  to  prescribe  large  temperature  changes. 

In  this  case  the  deformation-induced  temperature  changes  are  small  and  can 
be  neglected.  Then  R(t)  directly  determines  the  temperature  changes  and  we 
only  need  to  consider  Eq. (2). 

Simultaneous  changes  of  the  mechanical  variables  and  the  temperature 
occur  in  elevated  temperature  applications  and  are  used  in  mechanical  test¬ 
ing.  A  specific  example  is  thermal  fatigue  testing  [17-19] .  In  such  tests 
heating  is  accomplished  by  passing  an  electric  current  through  the  round  bar 
specimen.  Cooling  is  achieved  through  conduction  and  convection  [17-19],  and 
an  axial  (and  possibly  a  radial)  temperature  gradient  is  introduced  in  the 
specimen. 

Thermal  stresses  are  set  up  at  any  chosen  temperature  by  preventing 
the  axial  motion  of  the  specimen  using  a  suitable  clamping  device.  Clamping 
may  be  performed  at  the  minimum  or  maximum  temperatures  out  of  the  range  over 
Which  temperature  is  cycled,  and  different  stress-temperature  hysteresis  loops 
develop  depending  on  the  temperature  at  Which  clamping  takes  place  [17-19] . 
During  thermal  cycling  hardening  and/or  softening  can  occur  and  can  change  the 
hysteresis  loop  in  a  manner  similar  to  the  changes  observed  under  strain 
cycling  without  external  heating. 

We  idealize  this  test  situation  by  assuming  a  uniaxial,  homogeneous 
state  of  stress  and  adiabatic  conditions;  the  uniform  temperature  change  is 


prescribed. 


In  the  first  example  we  consider  a  bar  with  zero  initial  stress  cor¬ 
responding  to  an  initial  temperature  6q  and  constrained  to  have  zero  axial 
elongation  (e*0),  The  bar  is  subjected  to  adiabatic  heating  at  various 

constant  temperature  rates.  Figure  7  shows  the  compressive  stress  responses 

-2  +3 

predicted  at  heating  rates  from  10  to  10  degrees  K  per  second.  The  in¬ 
sert  of  Fig. 7  shows  the  temperature-dependent  elastic  modulus  E[0]  used  in 
this  simulation.  We  see  in  Fig. 7  that  the  stress  response  is  initially 
linear  elastic  in  the  temperature  change  but  subsequently  becomes  nonlinear 
and  "thermal  softening"  sets  in.  A  different  choice  for  the  temperature 
dependence  of  E[9]  will  alter  the  curvature  of  the  stress-temperature  curves 
and  the  thermal  softening  characteristics.  If  a  constant  modulus  is  chosen, 
the  rate  sensitivity  of  the  response  will  be  maintained  but  the  curves  will 
have  positive  slope  for  all  9  and  the  thermal  softening  represented  in  Fig. 7 
will  disappear. 

To  simulate  thermal  cycling  we  heat  to  either  900° K  or  1200° K  and  then 
cycle  between  900° K  and  1200° K  at  a  frequency  of  1  Hz*.  In  one  test  we  hold 
the  total  strain  constant  at  the  minimum  temperature  of  the  cycle  (900° K), 
Fig. 8.  Alternatively  in  Fig. 9  the  test  is  begun  at  the  maximum  cycle  temp¬ 
erature  of  1200° K.  In  each  of  these  cases  the  E[0]  of  Fig. 7  and  the  update 
parameters  of  Fig.l  are  used. 

We  see  that  in  each  case  stress-temperature  hysteresis  loops  develop 
similar  to  the  ones  observed  in  experiments,  see  Figs. 4  and  5  of  [19].  The 
increase  in  tensile  stress  shown  in  Figs. 8  and  9  is  higher  than  the  one 
observed  in  experiments.  This  increase  depends  on  the  updating  of  the 
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This  frequency  is  faster  than  the  one  used  in  experiments 


parameter  Xf  which  has  been  arbitrarily  selected,  flie  rounded  tips  of  the 
hysteresis  loops  are  due  to  the  imposed  sinusoidal  temperature  variation. 


Simultaneous  Mechanical  and  Thermal  Loading 

In  some  experiments  the  uniaxial  specimen  is  subjected  to  a  specified 

simultaneous  mechanical  and  thermal  loading  [20,21].  Here  we  consider  the 

case  of  a  uniaxial  deformation  field  corresponding  to  a  bar  which  is  cycled 

with  1  Hz  strain  control  (sine  wave)  while  being  rapidly  internally  heated. 

Here  internal  heating  corresponds  to  external  specification  of  the  internal 

heat  supply  R(t)  so  that  a  temperature  rate  G  =  300°K/sec  is  maintained 

throughout  the  strain  cycling  which  starts  at  an  initial  temperature  of 

9  =  300°  K.  Figure  10  shows  the  corresponding  stress-strain  response 

o 

throughout  this  combined  cycling  and  heating.  At  the  end  of  the  fourth 
cycle  the  temperature  of  1500°K  is  reached.  We  see  in  Fig. 10  that  this 
increase  in  temperature  distorts  the  hysteresis  loops  and  causes  cyclic 
softening  although  the  storage-updating  procedure  is  the  same  as  in  Fig. 1 
Where  cyclic  hardening  occurs.  The  distortion  and  softening  is  entirely 
due  to  the  large  specified  increase  in  temperature  and  the  temperature 
dependence  of  the  g- function  through  E[0].  Hie  temperature  dependence  of  E 
is  also  responsible  for  the  decrease  of  the  unloading  slopes  shown  in  Fig. 10. 

Discussion 

We  have  examined  this  theory  of  thermoviscoplasticity  in  the  case  of 
homogeneous  uniaxial  deformation.  Adiabatic  strain  cycling  and  also  temp¬ 
erature  cycling  have  been  simulated.  Also  studied  were  the  responses  in 
tensile  tests  at  different  strain  (stress)  rates,  the  responses  to  jump 
increase  and  jump  decrease  in  the  strain  (stress)  rate,  and  the  mechanical 
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response  to  heating  at  different  temperature  rates  (zero  elongation) .  For 
consistency  only  one  set  of  functions  k[  ]  and  g[  ]  have  been  used  through¬ 
out  these  simulations;  also  the  same  set  of  material  constants,  the  same 
value  of  E  when  taken  to  be  constant,  and  the  same  function  E[6]  when  taken 
to  be  temperature  dependent  have  been  used  throughout,  see  Table  1.  This 
procedure  facilitates  comparison  of  the  figures  and  assessment  of  the  theory. 

In  Figs.l-  5  we  see  that  the  temperature  changes  due  to  self  heating 
are  small  at  the  strain  (stress)  rates  used.  If  the  adiabatic  conditions 
were  to  be  removed  these  temperature  changes  could  easily  be  equilibrated 
through  heat  conduction. 

In  the  adiabatic  cycling  the  theory  presently  assumes  that  the  mechanical 
work  (the  area  of  the  hysteresis  loop  measured  from  zero  stress  to  zero  stress) 
is  converted  into  heat,  see  Eq. (4) .  Since  cyclic  hardening  and  softening 
involve  small  changes  in  the  area  of  the  hysteresis  loop,  hardening  and 
softening  do  not  have  a  significant  effect  on  the  temperature  change  gener¬ 
ated  per  cycle.  In  adiabatic  tensile  deformation  increasing  the  strain 
(stress)  rate  increases  the  steady-state  stress  (strain)  response;  this 
results  in  increased  self-heating  and  greater  temperature  change,  as  shown 
in  Figs. 4, 6.  In  cycling  increasing  the  frequency  increases  the  area  of  the 
hysteresis  loop  moderately  which  results  in  a  small  increase  in  the  tempera¬ 
ture  change  generated  per  cycle.  However,  in  a  fixed  period  of  time  a  sig¬ 
nificant  temperature  increase  would  result  from  increasing  the  frequency 
because  of  the  increased  number  of  cycles  completed  in  the  specified  time 
period. 

In  the  absence  of  a  temperature  change  the  mechanical  response  predicted 
by  Eq. (2)  is  symmetric  with  respect  to  the  origin.  In  tension  the  temperature 
initially  decreases  but  subsequently  increases,  Figs. 4, 6.  However,  in 
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compression  the  initial  temperature  response  is  immediate  and  continued 
heating,  see  Eq. (16)  here  and  Fig. 4  of  [8],  Thus  the  theory  predicts  that 
at  equal  strain  (stress)  magnitudes  the  temperature  differs  slightly  in 
the  respective  cases  of  tension  and  compression.  However  these  differences 
in  the  temperature  changes  are  small  (<  1°K)  and  result  in  negligible  dif¬ 
ferences  between  the  mechanical  responses  in  tension  and  compression. 

lhe  temperature  dependence  of  material  parameters  does  not  significantly 
influence  the  response  unless  large  temperature  changes  are  generated. 

These  can  occur  through  very  fast  adiabatic  tensile  (compressive)  loading, 
through  rapid  cycling,  or  through  internal  heating  (R  specified).  Tempera¬ 
ture  dependent  material  parameters  cure  important  in  thermal  loading  because 
of  the  large  temperature  increases  that  occur.  For  these  conditions 
significant  decreases  of  the  elastic  modulus  and  the  yield  strength  are 
observed  in  real  materials  which  must  be  modelled  by  the  theory. 

In  our  modelling  of  thermal  softening  the  only  temperature-dependent 
material  parameter  used  is  E[0]  which  is  shown  in  the  insert  of  Fig. 7.  Using 
this  E[6]  and  the  g-function  proposed  in  [10],  the  stress-strain  response  is 

predicted  to  have  an  initial  slope  which  significantly  decreases  with  increas¬ 

ing  initial  temperature.  Also  the  corresponding  steady-state  flow  stress  is 
decreased  considerably  as  the  initial  temperature  is  increased,  see  Fig. 8 
of  [8] .  The  temperature  dependence  of  the  modulus  is  responsible  for  the 
thermal  softening  present  in  the  thermal  loading  of  Fig. 7.  This  also  results 

in  the  distortion  of  the  strain  cycling  hysteresis  loop  of  Fig. 10.  Also 

the  decrease  in  the  flow  stress  upon  reloading  in  this  figure  is  due  to  the 
thermal  softening  which  overcomes  the  mechanical  hardening  present  in  Fig.l. 
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Figure  8  shows  a  decrease  in  the  flow  stress  on  the  first  compressive 
loading  which  is  much  more  pronounced  than  the  corresponding  case  in  Fig. 9. 
The  curvatures  in  the  tensile  and  compressive  segments  of  the  hysteresis 
loops  of  Fig. 9  are  different.  These  features  result  from  the  temperature 
dependence  of  E  and  the  thermal  softening  . 

We  have  simulated  adiabatic  strain  cycling  (Fig.l),  thermal  cycling 
under  constraint  (Figs. 8,  9),  and  mechanical  cycling  with  simultaneous  large 
internal  heating  (Fig. 10) .  The  same  g-function  augmentation  procedure 
(storage  and  updating)  and  the  same  successive  parameter  values  (Xf  from 
Table  1)  have  been  used  in  all  cases.  Although  parameter  changes  are  kept 
the  same  to  facilitate  comparison  between  the  figures,  the  strain  range  in 
Figs.l  and  10  is  1.8%  while  in  Figs. 8  and  9  the  range  of  the  mechanical 
strain  Z  is  approximately  .4%. 

In  real  materials  the  amount  of  hardening  will  decrease  with  decreasing 
strain  range;  smaller  increases  in  the  values  of  Xf  upon  updating  should  be 
used  in  Figs. 8  and  9  than  Figs.l  and  10  to  match  real  material  behavior. 

This  would  result  in  a  reduction  of  the  amount  by  which  the  tensile  stress 
increases  with  each  cycle.  Since  the  points  at  which  overstress  vanishes 
in  Figs. 8  and  9  correspond  to  different  temperatures  from  the  cases  of  Figs.l 
and  10,  a  temperature  dependence  of  the  parameter  Xf  in  addition  to  the  cycle 
dependence  of  Xf  would  introduce  additional  flexibility. 


* 

As  stated  earlier  the  tips  of  the  hysteresis  loops  are  round  due  to  the 
use  of  a  sine  wave  as  the  input;  they  would  be  sharp  in  the  case  of  a  saw¬ 
tooth  input. 


Hie  theory  presented  here  is  capable  of  representing  both  creep  and 
relaxation  behavior  corresponding  to  the  respective  test  conditions  a  =*  0 
and  e  *  0.  Hie  temperature  changes  predicted  during  such  deformations  are 
very  small  and  thermomechanical  coupling  is  not  important;  the  properties  of 
the  isothermal  theory  [7,9,23]  apply  in  this  case.  The  theory  predicts  that 
relaxation  terminates  at  a  nontrivial  equilibrium  value  of  stress  which  is 

determined  by  the  value  of  the  equilibrium  function  g[  ],  see  [7],  Primary 
and  secondary  creep  can  also  be  modeled  [7,23], 
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TABLE  1 


MATERIAL  FUNCTIONS  AND  PROPERTIES  USED  IN  THE  NUMERICAL  EXPERIMENTS 


1.  General  Parameters: 

a  =  14.4  X  lo"6  (°K)_1  pc  *  3.0  MPa/°K 

E  =  120  GPa. 


In  the  case  of  temperature-dependent  modulus,  E[0]  is  represented 
by  the  function  shown  in  Fig. 7. 


2.  Tlie  g[  ] -Function;  r,  Et,  Xf  are  material  constants,  see  [10] 

(Ete]-E  )  rtTi I 

g[Z,6]  -  EtZ  +  2r  tanh  [R  Xf  -  3]  loge  FcoshfV)] 


where 

U  =  R(X_  +  Z)  -  3  R=R  =3.6/X 

r  min  r 

V  =  R (Xf  -  Z )  -  3  Et  =  2.5GPa 

Xf  -  .003. 


In  cycling,  Figs.l  and  8-10,  the  parameter  X^  is  updated  at 

successive  occurrences  of  vanishing  overstress.  Successive 
values  of  Xf  are:  .00300,  .00600,  .00625,  .00650,  .00675,  .00700, 

0.00725,  .00750. 


3.  The  Function  k[  ]  is  Chosen  from  [24]  as: 

k-B  •  EXP[21.275  *EXP[-|a-g|/A]] 
B-  .2296  X  10~3  s 


FIGURE  CAPTIONS 


Figure  1 

Figure  2 
Figure  3 

Figure  4 

Figure  5 

Figure  6 
Figure  7 

Figure  8 

Figure  9 


Completely  Reversed  Adiabatic  Strain  Cycling  Using  a  1  Hz  Sine- 
Have.  Material  data  from  Table  1;  E  is  constant,  8  -  300° K. 

O 

Deformation- Induced  Temperature  Change  Corresponding  to  Fig.l. 

Tensile  Stress-Strain  Diagrams  for  Piecewise  Constant  Strain 

Rate  Loading  under  Adiabatic  Conditions.  The  equilibrium 

stress-strain  curve,  the  elastic  response  and  the  responses 

to  sudden  changes  in  strain  rate  are  also  shown.  Material 

data  from  Table  1;  E  =  const. ,  8  «*  300° K. 

7  o 

Adiabatic  Temperature  Response  for  the  Piecewise  Constant 
Strain  Rate  Tests  2-7  in  Fig. 3. 

Tensile  Stress-Strain  Diagrams  for  Piecewise  Constant  Stress 

Rate  Loading  under  Adiabatic  Conditions.  Hie  equilibrium 

stress-strain  curve,  the  elastic  response  and  the  responses 

to  sudden  changes  in  stress  rate  are  also  shown.  Material 

data  from  Table  1;  E  =  const,  8  *  300°  K. 

7  o 

Adiabatic  Temperature  Response  for  the  Piecewise  Constant 
Stress  Rate  Tests  2-7  in  Fig. 5. 

Adiabatic  Uniform  Heating  of  a  Constrained  Bar  (e  «  0)  at 
Various  Temperature  Rates.  Material  data  are  given  in 
Table  1.  The  E[8]  shown  in  the  insert  is  used;  8o=  300°K. 

Adiabatic  Thermal  Cycling  Between  900  and  1200°K  of  a  Bar 
Clamped  at  900°K  (e = const. )  Using  a  1  Hz  Sine  Wave.  E[8] 
of  Fig. 7  is  used  together  with  material  data  from  Table  1. 

Same  Temperature  Cycling,  Same  Material  Data,  and  Same  E[6] 

as  in  Fig. 8  Except  for  Clamping  which  Occurs  at  1200° K  (e = const. ) 

Simultaneous  1  Hz  Strain  Cycling  and  Uniform  Heating  at 

.  o  -1 

6  *  300  K  s  under  Adiabatic  Conditions.  Material  data 
from  Table  1  and  E[8]  from  Fig. 7  are  used. 
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Figure  10 


TEMPERATURE  CHANGE 
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Figure  5 
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